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We propose a scheme to generate temporal vector optical solitons in a lifetime broadened five- 
state atomic medium via electromagnetically induced transparency. We show that this scheme, 
which is fundamentally different from the passive one by using optical fibers, is capable of achieving 
distortion-free vector optical solitons with ultraslow propagating velocity under very weak drive 
conditions. We demonstrate both analytically and numerically that it is easy to realize Manakov 
temporal vector solitons by actively manipulating the dispersion and self- and cross-phase modula- 
tion effects of the system. 

PACS numbers: 42.65.Tg, 42.50.Gy 
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The vector nature of light propagating in a nonlin- 
ear medium has led to the discovery of a novel class 
of solitons, i. e. vector optical solitons, which are the 
solutions of two coupled nonlinear Schrodinger (NLS) 
equations describing the envelope evolution of two po- 
larization components of an electromagnetic field. In re- 
cent years, considerable attention has been paid to the 
temporal [E 0, i, i, i, Q and spatial 0, S, i, [Tfl] vector op- 
tical solitons in various nonlinear systems. Due to their 
remarkable property, vector optical soHtons have promis- 
ing appUcations for the design of new types of all-optical 
switches and logic gates [12]. 

Up to now, most vector optical solitons are produced 
in passive media such as optical fibers 0, i, i, i, 0, 0, i, 
H, , in which far-off resonance excitation schemes are 
employed in order to avoid unmanageable optical atten- 
uation and distortion. However, due to the lack of dis- 
tinctive energy levels, the nonlinear effect in such passive 
media is very weak, and hence to form vector solitons a 
very high input light-power is required. In addition, the 
lack of distinctive energy levels and transition selection 
rules also makes an active control very difficult. In par- 
ticular, it is hard to realize Manakov ^Hj temporal vector 
optical solitons in optical fibers because the ratio between 
self-phase modulation (SPM) and cross-phase modula- 
tion (CPM) is not unity and there is also detrimental 
energy exchange between two polarization components 
due to the existence of four-wave mixing effect. Manakov 
vector optical solitons are of great interest, not only be- 
cause the coupled NLS equations describing them have 
beautiful mathematical properties but also such solitons 
may be used to realize all-optical computingjl^. Differ- 
ent from spatial Manakov vector optical solitons, which 
have been observed more than ten years agoQ, temporal 
Manakov vector optical solitons have not been realized 
in experiment up to now. 

In this Letter, we propose a scheme to generate tem- 
poral vector optical solitons in a coherent five-level 
atomic system via electromagnetically induced trans- 
parency (EIT). This resonant EIT medium has been re- 
cently used to realize polarization qubit phase gate [14 
and reversible memory devices for photon-polarization 



qubit [15|. We show that two continuous- wave (CW) con- 
trol laser fields established prior to the injection of a 
pulsed probe field induce a quantum interference effect, 
which can suppress largely the absorption of the two 
orthogonal polarization components of the probe field. 
The scheme suggested here is fundamentally different 
from the passive ones due to the existence of distinc- 
tive energy-levels that make an active manipulation on 
the dispersion and nonlinear effects of the system pos- 
sible. In addition, contrary to all passive schemes the 
vector optical solitons produced in the present system 
may have ultraslow propagating velocity and their pro- 
duction needs only very weak input power. Furthermore, 
the controllability of the present scheme allows us also to 
realize easily temporal Manakov vector optical solitons 
by actively adjusting the parameters of the system. No- 
tice that scalar ultraslow optical solitons in EIT media 
have been investigated rece ntly[ll[I3, However, up 
to now there has been no study on the ultraslow vec- 
tor optical solitons in an active optical medium. Our 
study represents the first work in this direction and the 
results may have potential application in optical infor- 
mation processing and engineering. 

Consider a life-timed broadened five-level system (e. 
g. a Zeeman split atomic gas) interacting with a weak, 
pulsed, linear-polarized probe field of central frequency 
lOp/{2tt) and two strong, linear-polarized CW control 
fields of frequencies a;ci/(27r) and a;c2/(27r), respectively. 
The two polarization components of the probe field drive 
respectively the transitions from |3) ^ |2) and |3) ^ |4), 
while the two control fields drive respectively the tran- 
sitions from |1) ^ |2) and |5) ^ |4) (see Fig. 1(a)). 
The atoms are trapped in a cell with the temperature 
lowed to 0.5 /iK to cancel Doppler broadening and colli- 
sions. A possible arrangement of experimental apparatus 
is shown in Fig. 1(b). The electric-field of the system 
can be written as E=(e+fp+ -|-e_fp_) exp[i(fcpZ — Wpi)] -I- 
Gcifci eyiY>[i{kciy-uJcit)]+ec2£c2 exp[i{kc2y-^c2t)]+c.c.. 
Here e+ = {'k+iy) / \/2 (e_ = (x— iy)/\/2) is the probe-field 
unit vector of the ) circular polarization compo- 

nent with the envelope £p+ {£p- ) , which drives the tran- 
sition |2) ^ |3) (|3) ^ |4)). Gel (ec2) is the unit vector 
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FIG. 1: (a): Energy level diagram and excitation scheme of a 
five-level atomic system interacting with a weak, pulsed probe field 
of Rabi frequency Qp and two strong, CW coupling fields of Rabi 
frequency Qcl and Qc2, respectively, (b): Possible arrangement of 
experimental apparatus. 



FIG. 2: Absorption spectra of flpi (solid line) and Qp2 (dashed 
Une) with parameters r2 ~ r4 = 27r X 6 MHz, Ti ~ Ta ~ Ts = 2 
kHz, K32 ~ /t34 = 1.0 X 10'' cm-is-\ Ud = Qc2 = 2.0 x 10^ 
s~\ and (5ci = 5^2 = 0, and A = 2.0 X 10*^ s'^ . The dash-dotted 
(dotted) line represents the absorption spectra of Qpi {Qp2) when 
the control fields are switched off. 



of the control field with the amplitude Ed (Ec2), which 
drives the transition |1) ^ |2) (|4) <^ |5)). Thus the 
system is composed of two EIT A-configurations, both of 
them share the ground-state level |3)[l4lT5j. 

In interaction picture, the atomic response of the sys- 
tem under rotating-wave approximation is described by 



d_ 

dt 



idi Ai = -in*^A2 



(la) 



— + id2n ) A2n = -if^cn ^4n-3 " «f^pn^3, (lb) 



(Ic) 



— + ids ] A5 = -ifll^A^, 



{n—1,2) with J2j=i l^jP— where Aj is the probabil- 
ity amplitude of the bare atomic state \j) (with eigenen- 
ergy Cj), f7pi=-(p23-e_)£'p_/fi, Op2=-(p43-e+)fp+/ft, 
^ci=-{P2i-^ci)£ci/h and r^c2=-(P45-ec2)fc2/fi- are half 
Rabi frequencies with p^^- being the electric dipole ma- 
trix element associated with the transition from \ j) and 

In Eq. ^ we have defined di = {5p — Sci) — 
d2=Sp- ir2/2, d3 = -ir3/2, d4 - {Sp + A) - ir4/2 and 
d^ = (dp + A - Sc2) - jr5/2 with Sp (£2 - £3)/^ - Wp, 
Sci = (e2 - ei)/h - Wci, and Sc2 = (£4 - £5)/^ - t^c2- 
Tj is the decay rate of the state A — {2^B/K)gB is 
the Zeeman shift of the upper atomic sublevel with /i^ 
the Bohr magneton, g the gyromagnetic factor and B the 
applied magnetic field. 

The equation of motion for r2p„(z,t) can be obtained 
by Maxwell equation under slowly-varying envelope ap- 
proximation 



d_ 1 5_ 
dz ^ cdt^ 



K3Z^2„^3 = 0, (/ = 2n) (2) 



where K32 = A/'a|p32 • e-\'^ujpl{2heQc) and K34 = A/'q|p34 • 
£+pa;p/(2?i£oc) with A/'a being the atomic density, £0 the 
vacuum dielectric constant and c the light speed in vac- 
uum. 

Before solving Eqs. ([!]) and ([2]), we examine the lin- 
ear properties of the system, which provide main con- 
tributors to pulsed spreading and attenuation. We as- 
sume that the probe field is weak so that the atomic 



ground state |3) is not depleted, i.e., A3 wl. Taking 
U,pn and Aj [j = 1,2,4,5) as being proportional to 
exp[i(fc(a;)z — Lot)]^ one can get two branches of linear 
dispersion relation ki(uj)=uj / c -\- K32(i^ — di) / Di{uj) and 
k2{uj)=uj/c + K34(w — d^) / D2{oj), corresponding to a" 
and (7+ components of the probe field, respectively. Here 
we have defined £'i(a')=|riciP — {io — di)(w — ^2) and 
£)2(ti^)=|r2c2p — [oJ — di){LO — d^). Shown in Fig. 2(a) is 
the absorption spectra of Q.pi and Q.p2 . We see that near 
the central frequency of the probe field (i. e. at u}=Q), 
both Im(fci) and Im(fc2) are close to zero, which means 
the absorption is greatly suppressed. The reason of such 
suppression of the probe field is due to the introduction 
of the two strong control fields that induce a quantum 
interference effect and thus make the two polarization 
components of the probe field transparent. 

Because of the existence of dispersion effect, the probe 
field will distort during propagation. In the following we 
show that the SPM and CPM effects of the system may 
balance the dispersion. To this aim we apply the method 
of multiple-scales [TtI [isj to investigate the weak nonlin- 
ear evolution of the probe field. We make the asymptotic 

expansion Aj =X]l^o A^'^j''' ^'^'^ f^pn=X)i^i m'^p" with 
A^^^ — \ and A^^—Q (^=1,2,4,5), where /i is a small pa- 
rameter characterizing the small population depletion of 
the ground state and all quantities on the right hand side 
of the asymptotic expansion are considered as functions 
of the multi-scale variables z/=/x'z and ti=^H. 

The leading order solution is given by fl^n = 
Fn exp{z[fc„(a;)zo — wio)]} with F„ being yet to be de- 
termined envelope functions. At the second order, a 
solvability requires i[dFn/dti + VgndFn/dzi] — 0, where 
l/g„=l/i^i„ with Kn = l/c+«32(|f^ciP + rf?)/(|r!cip- 
did2 f and K12 = l/c+K34(|r!c2p + di)/(|r!c2p-rf4rf5)'. 
The solvability condition at the third order yields two 
coupled NLS equations 



K2nd^Fn 



dz2 



dtl 



-(W„„|F„|2-fW'„„,|F,„ne-2"'.-^F„ = 0, 

(3) 
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(m,n=l, 2] m ^ n) where 



K.2 



-2^3; 



(|d4„'-3|' + |^^c„'|') 



(n, n'=l, 2) are the coefficients characterizing the dis- 
persion {K2n), SPM (Wnn) and CPM {Wnm, u ^ m) of 

the two polarization components of the probe field, re- 
spectively, an = ^J.'^an with a = Im[Xo„] [ifo„=A:„(0)]. 
When returning to original variables and introducing 

S=il/Vg, - l/l^g2)/2, Vg = 2VglVg2/{Vgl + ¥^2) , aUd T = 

t — z/Vg, Eq. ([31) can be written as the dimensionless 
form 



d 



ds 

{^Qnn I 



+ 9An Un + (-1)' 



dUn 

97 



2 9ct2 



0, 



(4) 



where s=z/L]:), ct—t/tq, and u„=(f2p„/[/o)e 
(Ko„=Re[ii'on]), gAn^ctnLo, gs^sign{S)LD/Ls, 5di = 
-^■21/1^^22!, 5D2=sign(/4:22) and 5(nn'=W„«'/|W^22|- Here 
we have defined Li:)=Tq /\K22\ (dispersion length), and 
^i5=''o/|'^| (group velocity mismatch length), tq is typi- 
cal pulse length of the probe field. In order to get soliton 
solutions we have assumed that typical nonlinear length 
Lnl — 1/(C^o|W^22|) is equal to the dispersion length Ld. 

Because coupled NLS Eq. ([4]) has complex coeffi- 
cients, generally a vector soliton does not exist. How- 
ever, as we show below that practical parameters can 
be found based on the EIT effect and hence the imagi- 
nary part of the coefficients can be much smaller than 
the corresponding real part. This leads to a shape- 
preserving vector optical soliton solution that can prop- 
agate for an extended distance without significant de- 
formation. The system admits bright-bright, bright- 
dark, and dark-dark vector soliton solutions through a 
balance between the dispersion and nonlinear effects. 
The bright-bright vector soliton solution reads w„ = 
V„sechcrexp[i(7^nO' + Qns)] (n = 1,2) if the parameters 
fulfill the condition 922901 = 9i29D2- Here we have de- 
fined Vn = {-l)''-^9s/9Dn, Qn = -9y{29Dn) " 9Dn/2, 

and V2 = [(.9m — 5ii"^'i)/5i2]^^^- A bright-dark vector 
soliton solution is given by ui = Visechcr exp [1(7^1 cr + 
Qis)], U2 = V2 tanh(Texp[«(P2CT + Q2S)], where Vn = 
{-lY-^gs/9Dn, Qi = -Vigs - - V!)/2 - .912 V|, 

Q2 - V2gs + 9D2Vi/2 - g22Vl and V2 = [{gnVf - 
9d\)I gi-i^^"^ ■ Here Vi is a free parameter. 

Now we show that a realistic atomic system can be 
found that allows the bright-bright vector optical soliton 
described above. We consider a cold alkali atomic vapor 
(e. g., rubidium or cesium atoms) with the decay rates 
r2 ~ r4 = 0.5 X 10^ s-\ and Ti ~ Tg ~ Eg = 1.0 x 10^ 
s-i. We take K32 ~ K34 = 1-0 x 10^ cm-^s^i {Ma 10^" 



cm 



^), Vl^x = 0.^2 = 1.6 X 10** s-i, (5p = 1.0 X 10** 



A = 2.0 x 10" s-\ = 0, and 6^2 
With the above parameters, we obtain Kq\ 



3.0 X 10^ 



s 

-6.41 



O.lOi cm-i, Kq2 = -6.38 -t- O.lOi cm"!, ifn = (14.62 - 
0.47i) X 10"** cm-^s, Kx2 = (14.72-0.47i) x 10"** cm-^s, 
K2X = (-4.56 + 0.25i) X lO^^^ cm-^s^, K22 = (-4.64 -|- 
0.26i) X 10 



-15 



cm 



= (-9.37 + 0.15i) X 10- 
Q-ls^ 1^12 = (-9.44-H0.15i) X 10-16 cm-^s^, W2X 
-9.34 -h 0.15«) X 10-16 cm-^s^, and 1^22 = (-9.40 
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Notice that the imaginary parts 



(■ 

0.15i) X 10-16 cm 
of these quantities are much smaller than their relevant 
real parts. The physical reason for such small imaginary 
parts is due to quantum destructive interference induced 
by two CW control fields (i. e. EIT effect). We obtain 
Ls = 116.8 cm and = 0.8 cm with tq = 6.0 x 10"® 
s and J7o = 3.7 x 10^ s^^. The dimensionless coefficients 
read 91, = 0.007, gui = —0.98, gD2 = —1.0, and 511 ~ 
gi2 ~ 521 — 322 = —1.0. The group velocities of the 
two polarization components are respectively given by 
Re(Vgi) = 2.28 X 10-4 c and Re(yg2) = 2.26 x 10-^ 
c, which means that the two polarization components of 
the vector optical soliton propagate with nearly matched, 
ultraslow propagating velocities. 

As we have stressed, different from the passive media 
such as optical fibers @, i, 0, i, i, 0, i, M, El the pa- 
rameters of our present EIT medium can be actively ma- 
nipulated. Consequently the coefficients of Eq. (|4]) can 
be easily adjusted to allow us to realize a Manakov sys- 
tem, which is a completely integrable and can be solved 
by inverse-scattering transform[ll]. In fact, with the pa- 
rameters given above Eq. ([4]) can be written as a near 
Makakov system {m,n=l, 2; m ^ n) 



ds 



1 d'^Un 



2 



+ (k«P + l"mP)Wn = Rn{Un), (5) 



with i?„(u„) ~ — 0.08m„ describing the linear absorp- 
tion effect. We see that i?„ is indeed a small quantity 
which can be taken as perturbation. The vector soli- 
ton solution of Eq. ([5]) after neglecting i?„ is ui — 
cos0sech(cT) exp(is/2) and U2 = sin0sech(a) exp(is/2) 
with being a free parameter. Note that since the in- 
jected probe field is tt (i. e. linear) polarized, the two 
polarization components should have equal amplitude, i. 
e. 9 = 7r/4. 

Shown in Fig. 3 is the evolution of the two polariza- 
tion components of the probe field versus dimensionless 
time i/ro and distance z/ Ljj. The plots are obtained by 
numerically integrating Eq. ^ by using a split-step fast 
Fourier transform method and the bright-bright soliton 
solution given above as an initial condition. To demon- 
strate the balance between the dispersion and nonlinear 
effects, we change the probe field amplitude Uq while 
keep other parameters the same as those given above. 
Fig. 3(a) shows the case when the dispersion is domi- 
nant over nonlinearity (i. e. UqTq < \/\K22/W22\). We 
see that in this case both polarization components spread 
and attenuate seriously. Fig. 3(b) shows that case when 
UoTo—y^\K22/ W22 1 , i- e. there is a balance between the 
dispersion and nonlinearity. In this situation a shape- 
preserving propagation of vector optical soliton over long 
distance (z=3.2 cm) is achieved. 
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FIG. 3: Evolution plots of the two polarization components of the 
probe field, (a): Dispersion dominant case with Uo=3.7 X 10® s~^. 
(b): The case of a balance between dispersion and nonlinearity with 
Uo=3.7 X lO'^ s~^. Brighter shading marks higher Rabi frequency. 
The propagation distance is z = ALo = 3.2 cm. 




FIG. 4: The bright-bright vector soliton evolution obtained by 
integrating directly from Eq. JTJ and JSJ without any approxima- 
tions. The solid (dashed) line is for the a— {cr+ ) component of the 
relative probe intensity \Qpi/Uq\'^ {\^p2/Uo\^)- The propagation 
distance is z = l.SLo = 1.2 cm. 

The input power of the vector optical sohton can 
be calculated by Poynting's vector. It is easy to get 
the average flux of energy over carrier-wave period 
P„=P„j-„ax sech^ [(t — z/Vgn)/TQ], with the peak power 



Amax ~ Amax = 8.9 X lO""* HiW. Here we have taken 
IP23I ~ IP43I = 2.1 X 10~^^ cm C and the beam radius of 
the probe laser R± = 0.01 cm. We see that to generate an 
ultraslow vector optical soliton in this active system only 
very low input power is needed. This is drastically differ- 
ent from the vector optical soliton generation schemes in 
passive midia[i, S, i, B H 0, B B Il3| where much higher 
input power is needed in order to bring out the nonlinear 
effect required for the soliton formation. 

To make a further confirmation on the vector soli- 
ton solutions obtained and check their stability, we have 
made additional numerical simulation directly from Eq. 
([1]) and ([2]) without using any approximations. Fig. 4 
shows the simulation result by taking m„ (z = 0) = 
2sech(2-\/2cr) (n = 1, 2) as initial condition. We see that 
the soliton radiates a small part of energy in its tail but 
is fairly stable during propagation. 

In conclusion, we have proposed a scheme to create 
temporal vector optical solitons in a coherent five-level 
atomic system. Such solitons can have ultraslow propa- 
gating velocity and may be produced with extremely low 
input power. We have demonstrated both analytically 
and numerically that it is easy to realize Manakov tem- 
poral vector optical solitons by actively manipulating the 
dispersion and nonlinear effects of the system. Due to the 
robust propagation nature, the ultraslow vector optical 
solitons suggested here may have potential application 
in optical information processing and engineering under 
a weak-light level. 
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